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The dynamics of a point charged particle which is driven by a uniform external electric field 
and moves in a medium of elastic scatterers is investigated. Using rudimentary approaches, we 
reproduce, in one dimension, the known results that the typical speed grows with time as t^^^ and 
that the leading behavior of the velocity distribution is e"'"' In spatial dimension d > 1, we 
develop an effective medium theory which provides a simple and comprehensive description for the 
motion of a test particle. This approach predicts that the typical speed grows as t^'^ for all d, 
while the speed distribution is given by the scaling form P(u,t) = {u)~^ f{u/{u)), where u = li'l^'^^, 
{u) ~ and f{z) oc z'''~^-'''^e~^ ^'^ . For a periodic Lorentz gas with an infinite horizon, e. g., for 
a hypercubic lattice of scatters, a logarithmic correction to the effective medium result is predicted; 
in particular, the typical speed grows as (t Int)^''"^. 

PACS Numbers; 02.50.-r, 05.40. +j, 05.60+w 
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I. INTRODUCTION 

At the turn of the century, Drude developed a quali- 
tative theory for electrical conduction in metals 0|. To 
establish a more solid basis for the Drude theory, Lorentz 
[H suggested an ideaUzed model for electron transport in 
metals in which: (i) electron-electron interactions are ig- 
nored, (ii) the background atoms are considered to be im- 
mobile spherical scatters, and (iii) the electron-atom in- 
teraction is described by elastic scattering. This Lorentz 
gas j^] has played a large role in developing our under- 
standing of diffusive transport in random media. 

An important feature of the Lorentz gas is the indepen- 
dence of the electrons. This implies that the underlying 
Boltzmann equation for the evolution of the electron ve- 
locity distribution function (VDF) is linear. Because of 
this simplification, the Boltzmann equation has proven 
fruitful in understanding the properties of the Lorentz 
gas (see, e. g., Q and references therein). These inves- 
tigations have established that, under relatively general 
conditions, a test particle moves diffusively and that its 
diffusivity can be computed in terms of the geometric 
properties of the background scatterers. The Lorentz 
model is also simple enough to be amenable to rigor- 
ous analytical studies (see, e. g., |]5|-|l0|). In particular, 
for a periodic Lorentz gas in two dimensions with an "in- 
finite" horizon (z. e., there exist free trajectories of infi- 
nite length), anomalous diffusion of the form (r^) oc t Int 
has been proved ; this anomalous diffusion is also ex- 
pected to hold in arbitrary dimension. 

Paradoxically, much less is known about the prob- 
lem which originally motivated the introduction of the 
Lorentz gas model, i. e., the motion of a charged test 
particle in a scattering medium under the influence of 
a spatially uniform electric field. Lorentz himself con- 
structed an approximate stationary solution to the Boltz- 
mann equation by a perturbative expansion around the 



Maxwell-Boltzmann distribution | 
Lorentz reproduced the basic resu 



2]. From this solution, 
ts of the Drude theory. 



Unfortunately, the starting point of Lorentz's analysis is 
erroneous. If the scattering is elastic (no dissipation), 
then an electron will necessarily gain kinetic energy as it 
accelerates in the field and a stationary asymptotic VDF 
will not exist. This dilemma motivated investigations of 
the field-driven Lorentz gas in which some form of dis- 
sipation is explicitly incorporated pd|-p^, so that it is 
possible to obtain Ohm's law. 

In the absence of dissipation, however, Piasecki and 
Wajnryb p4| were apparently the first to recognize 
the fundamental ramifications that arise from the non- 
stationarity of the system. From an exact solution to the 
Boltzmann equation in one dimension and an asymptotic 
solution for general d and with the crucial assumption 
of collision isotropy, they found: (i) the typical veloc- 
ity, frms, grows with time as t^/'^, and (ii) the VDF has 
a non-stationary, but symmetric asymptotic form whose 
controlling factor is e"'"' 

Our goal in the paper is to develop simple and physi- 
cally transparent approaches to understand the behavior 
of the field-driven Lorentz gas. We begin by considering 
the one-dimensional system in Sec. II, where we substan- 
tially reproduce the results of Piasecki and Wajnryb Q. 
We first construct a random walk argument to explain the 
mechanism that gives rise to the slow t^/^ increase of the 
root-mean-square velocity (wrms) with time. This argu- 
ment relies on the assumption that each scattering event 
is spatially isotropic. This isotropy is a basic feature of 
elastic scattering from an immobile sphere only for the 
physical case of three dimensions jl^]. If one postulate 
isotropic scattering even in one dimension, the results ob- 
tained can be expected to mimic three three-dimensional 
behavior. The Langevin and underlying Fokker-Planck 
equations for the speed distribution function (SDF) are 
also investigated to obtain both the time dependence of 
the typical speed and, more generally, the asymptotic 
form of the SDF. Finally, we develop a Lifshitz argument 
p6[ to reproduce the asymptotic form of the SDF with 
minimal calculation. 
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In Sec. Ill, wc study the field-driven Lorentz gas for 
arbitrary spatial dimension d. We argue that there are 
basic differences in the field-driven Lorentz gas for d = 1 
and for d > 1. In greater than one dimension, the 
freely-accelerated trajectory segments between scatter- 
ing events are biased by the field, leading to anisotropy 
in the spatial position of the test particle at the next 
scattering event. To account for this bias in a relatively 
simple manner, we introduce an effective medium the- 
ory. In this description, a particle begins at the center of 
a "transparency" sphere of radius equal to the mean- free 
path. The particle freely accelerates until it reaches the 
sphere surface. This defines a collision, whereupon the 
test particle starts at the center of a new transparency 
sphere. We will generally assume isotropic scattering in 
each collision, a feature of elastic scattering from hard 
spheres in three dimensions. However, there is preferen- 
tial backscattering for For d <3 and preferential forward 
scattering for d > 3. This short-range correlation appar- 
ently docs not affect the asymptotic motion of the test 
particle, so that we generally focus on isotropic scatter- 
ing. 

For isotropic scattering, it is simple to quantify the 
field-induced bias of the test particle as it moves within 
a transparency sphere. A random- walk argument of a 
similar spirit to that given in one dimension indicates 
that the influence of the bias is of the same order as the 
stochasticity caused by scattering. This implies that the 
SDF will obey one-parameter scaling. From the solu- 
tion to the underlying Fokker-Planck equation, we find, 
P{u,t) oc t-i/2z(''-i)/3exp(-zV2), where u = is 
a convenient "speed" variable and the scaling variable z 
is proportional to ujl}!'^. We also extend the effective 
medium theory to the case where the mean-free path 
is chosen from a distribution cx This form 

accounts for the asymptotic distribution of mean free 
paths encountered by a test particle when (small) scat- 
tering centers are on regular lattice sites. As might be 
expected, when /x > 3, corresponding to a finite second 
moment of the transport of the test particle is 
nearly identical to the case where the mean-free path is 
fixed. However for /i < 3, i. e., for a distribution p(t^ 
with = oo, faster asymptotic transport arises. The 
borderline case of = 3 corresponds to a lattice array of 
scatterers such that an infinite horizon arises, and loga- 
rithmic corrections in the transport laws arc predicted to 
occur. 

In Sec. IV, we present Monte Carlo simulation results 
for the motion of a test particle in the two-dimensional 
effective nicidium. When the radius of the transparency 
circle is fixed, we obtain excellent agreement between 
simulation results and our theoretical predictions for the 
case of isotropic scattering. For simulations based on the 
correct scattering law for hard circles in two dimensions, 
virtually identical results are obtained, i e., short range 
antipersistence appears to be asymptotically irrelevant. 
We also consider the case of a power-law distribution of 
radii for the transparency circle, p[£) oc While the 



physically interesting case of /i = 3 corresponds to the 
borderline of applicability of our naive effective medium 
approach, numerical results indicate transport properties 
which are close to those obtained for a fixed radius trans- 
parency circle. 

In Sec. V, we present a brief discussion and summary. 



II. LORENTZ GAS IN ONE DIMENSION 
A. Random Walk Argument for the RMS Velocity 

Consider a sizeless test particle (electron) which moves 
with constant acceleration a = eE/m, where e and m 
arc the electron charge and mass, and E is the electric 
field. The electron moves in a medium of equally spaced 
point scatterers with separation i. To mimic the behavior 
of a three-dimensional system with isotropic scattering, 
we allow the particle to hop with equal probability to 
its nearest neighbor on the left or the right after each 
collision. Thus the electron trajectory consists of freely 
accelerating segments which are punctuated by isotropic 
scattering events. This is simply an isotropic random 
walk, but with position- and direction-dependent time 
increments between successive steps. 

We use this physical picture to compute the behav- 
ior of the typical velocity as a function of time. Energy 
conservation gives 



-mv^ — eExn = const. 



(1) 



where Vn and x„ refer to the electron velocity and po- 
sition immediately after the n^^ scattering event. We 
rewrite this as 



2 2ei;, , , 2eEi 

-v^ = —{xn+i-Xn)=±—^. (2) 



Because of the postulated isotropic scattering, v'^j^-^ — 
is equally likely to be positive or negative. Thus we 
conclude that undergoes a simple random walk as 
a function of n, with an elementary step size given by 
Vq = 2eE£/m. As a result, 



or 



(3) 



To determine the dependence of Wrms on time, we write 
the time increment between successive collisions as 



(4) 



The last approximation applies when the typical speed is 
large so that the acceleration between scatterings can be 
neglected, an assumption which can be verified a poste- 
riori. The total elapsed time for n collisions is therefore 



dtk ~ — 



fe=i 



dk 



(5) 
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Solvine 
Eq. (i 



for n as a function of time and substituting into 
gives the fundamental result 



1/3 



Vot 

£ 



1/3 



(6) 



It is instructive to compare the time dependence of 
v-n-ns with that of the average velocity in the field direc- 
tion. The latter can be computed from the recursion 
relation 



±Vn H , 

Vn 



(7) 



By isotropy, the factor of ±1 occurs equiprobably for each 
scattering. Since the typical speed grows indefinitely, we 
again ignore the acceleration during the free flight be- 
tween adjacent sites, so that Vdritt = (vn) ~ {(ai/vrms))- 
As a function of time, this may be rewritten as 



Vdrih{t) ( — 



1/3 



(8) 



Thus the average drift velocity decreases with time, even 
though the rms velocity grows with time. Therefore the 
VDF becomes systematically more isotropic in the long 
time Hmit Q. 

Finally, making use of Eq. (^), one can estimate the 
average displacement {x{t)) in the field direction to be. 



{x{t)) ^Vd,i{t{t)tr^ {alhy/^ 



(9) 



Alternatively, this same result follows directly from en- 
ergy conservation, Eq. (|l|), and the time dependence of 
Vrms(i) from Eq. B). 



B. The Speed Distribution 

We now derive the speed distribution function using 
simple approaches which obviate the need to solve the 
Boltzmann equation, as given in jl4j. First consider the 
Langevin equation to describe how the typical speed de- 
pends on n. Since is randomly incremented or decre- 
mented by a fixed amount Vq in a single collision, wc may 



write, in the large-rt limit, 

<K 

dn 



vlT]{n), 



(10) 



where the noise has zero mean, {r]{n)) — 0, and is tem- 
porally uncorrelated, {'ri{n)rj{n')) = 5{n — n'). Since we 
are interested in the n — s- oo limit, the continuum re- 
sult for the above correlation function is appropriate. In 
this limit, the amplitude distribution of the noise is also 
Gaussian. Consequently, the Langevin equation yields a 
Gaussian distribution for with a dispersion equal to 



To determine the time dependence of the speed dis- 
tribution, we transform from n to t by writing dt = 



tdn/\vn\, so that 



dvl 
dn 



2(. 



d\Vn 

dt 



(11) 



Next, we transform the dependence of the noise correla- 
tion from n to t. Writing 5{n — n') = 5[t — t') gives 

{V{n)v{n')) = {r]{t)r]{t'))£/\v\, so that rj{n) = r?(t)VW- 
Substituting this into the Langevin equation, Eq. ([lC|), 
gives 




dt 



dt - aVI^^'^' 



(12) 



(13) 



Thus we conclude that the speed distribution function, 
P{u, t), is Gaussian in m = with a dispersion which 

is proportional to Vq/£. Then the VDF is determined 
from the identity P{v, t)dv = P{u, t)du to yield 



P{v,t) 



exp 



9eaH 



(14) 



An independent and appealing approach to obtain the 
VDF is by a Lifshitz tail argument |l^. This method 
is based on matching the assumed scaling form of the 
VDF with the "extreme" contribution that arises from a 
particle which is scattered in the field direction at each 
collision. This extreme tail can usually be estimated by 
elementary means, and matching this to the scaling form 
determines the VDF. Although this approach is heuristic, 
its advantages are simplicity and wide applicability. 

Our starting point is to assume that the VDF can be 
written in the scaling form 



Piv,t) 



1 



-/(w/Wrms), 



(15) 



where the scaling function f(z) is expected to approach a 
constant as z — > 0, and vanish faster than any power law 
for z — > oo. Generally, this large-z asymptotic behavior 
has the quasi-exponential form 



/(^) 



exp ( — z 



(16) 



which defines the "shape" exponent S of the distribution. 

Let us now consider the trajectory in which the test 
particle is perpetually scattered parallel to the field, so 
that its speed is simply v — at. Substituting into Eq. (|l^ ) 
and using Eq. (Bh for Vrms gives. 



P{v = at,t) 



(17) 
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On the other hand, the probabihty P„ that n scattering 
events are all parallel to the field equals 2~". For this uni- 
formly accelerated motion, the correspondence between 
n and the time is simply given by = ni. Thus 

writing P„ as a function of time and matching with the 
argument of the exponential in Eg. (|l7[) , gives (5 3, in 
agreement with the exact solution |14|]. Parenthetically, 
if we define a size exponent v through fims ^ , then the 
general scaling relation iQ between the size and shape 
exponents, S = (1 — fails for the field-driven Lorentz 
gas. 



III. LORENTZ GAS IN GREATER THAN ONE 
DIMENSION 

A. Effective Medium Approximation 

The field-driven Lorentz gas in greater than one dimen- 
sion presents a variety of theoretical and computational 
challenges. Numerical simulations of the dissipationless 
system are prone to large fluctuations and quantitative 
conclusions are not readily obtained Because of 

this computational difficulty and also because dissipa- 
tion arises in any physical realization of the Lorentz gas, 
simulation work has primarily focused on the field-driven 
system with dissipation. This is achieved by either allow- 
ing for inelasticity in collision events |]l2| , or by introduc- 
ing a "thermostat" which continuously extracts energy 
from the particle during its free motion to maintain a 
constant kinetic energy While much is known 

about these dissipative systems our interest is in 

the nonstationary behavior of dissipationless system. In 
particular, we wish to understand how the typical speed 
of a test particle grows with time and the form of the 
resulting SDF. 




X 



FIG. 1. "Transparency" sphere that surrounds a scatterer. 
After a scattering event, the test particle moves freely on 
a parabohc trajectory until the next collision at the sphere 
boundary. The initial and final angular position of the test 
particle, 8 and a, respectively, are indicated. The critical tra- 
jectory is defined by the condition that the final longitudinal 
position of the test particle is at a: = 0. 

Because of the inherent difficulties in describing the 
motion of a test particle in a regular lattice of scatterers, 
we introduce an effective medium approximation in which 
the true trajectory is replaced by an effective, but phys- 
ically equivalent, trajectory whose properties are readily 
calculable (Fig. 1). We assume that immediately after 
each scattering event, the test particle starts at the center 
of a transparency sphere of radius equal to the mean-free 
path £. The test particle then freely accelerates until the 
next collision when the surface of this sphere is reached. 
The collision point defines the center of the next trans- 
parency sphere. This construction is repeated to gen- 
erate a particle trajectory which consists of parabolic 
segments (the biased free-particle motion between col- 
lisions), which are punctuated by collision events. We 
assume isotropic scattering so that the outgoing particle 
direction is randomized. This isotropy actually occurs 
for hard-sphere scattering only in three dimensions. An 
elementary computation shows that there is preferential 
back scattering for d < 3 (with complete back scattering 
for d = 1), and preferential forward scattering for d > 3. 
However, this persistence for d > 3 or antipersistence 
for d < 3 appears to be asymptotically negligible (see 
below) . 



B. The Typical Speed 

To estimate the typical speed, we need to quantify the 
deflection of a trajectory during free flight. As we shall 
show, this leads to an effective bias which vanishes as 
the inverse square of the particle speed. Let us define 
trajectories whose collision points are in the hemisphere 
X > as positively biased and vice versa. Separating 
these trajectories is a "critical" trajectory, in which the 
next collision point is also at x = (Fig. 1). (This 
critical trajectory exists only if the initial speed satis- 
fies V > vo/V2] otherwise all trajectories are deflected 
towards increasing x) However, since the typical speed 
grows as a power law in time, the role of trajectories in 
which the speed at some stage is too small to define a 
critical trajectory is expected to be negligible. 

By elementary mechanics, the inclination angle of this 
critical trajectory is given by 

1 . _i / TT /uo\^ vo 
0=-sm T ~ TT + t;- , as ^ 0. (18) 

2 \ J 2 \2vJ V ^ ' 

The longitudinal motion of the test particle can thus be 
viewed as a biased one-dimensional random walk, with 
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the bias at each step proportional to e = (i'o/2w)^. Fol- 
lowing the same steps as those given in Eqs. (1) - (6), 
the velocity increment between scatterings is given by 
v'^_^_l — v^ — ±Wq, but with the ± sign now occurring with 
respective probabilities ^(l±ae), where a is a dimension- 
dependent number of order unity. Thus, in addition to 
the stochastic particle motion given in Eq. (|^), a deter- 
ministic contribution also arises. This latter component 
gives, for the n dependence of the speed, 



(19) 



To determine the time dependence, we relate n at t by 
using 



E" „ f"' dk i pn 



(20) 



Solving for n as a function of t, substituting into Eq. (19), 
and eliminating the factor e, we find that the speed is 
given by 



1/3 



(21) 



This is the same time dependence as one dimension, 
where there is no deterministic bias in the motion. For- 
tuitously, the manifestations of the isotropic scattering 
and the field bias are of the same order in our effective 
medium theory. This coincidence leads to a distribution 
of speeds which can be described by single parameter 
scaling. 



C. The Speed Distribution 

To determine the speed distribution, we adapt our 
approach used in one dimension. First, we derive the 
Langevin equation for the dependence of the typical 
speed on n, from which the underlying Fokker- Planck 
equation may be written and then solved. 

In a time At after collision n (and before collision 
n + 1), the particle will be at 



2a(r • e) = 2avT{n ■ e) -|- (ar)^ 



(24) 



Combining Eqs. ( |23| ) and (|24|), the time increment can 
be eliminated to give the analog of Eq. (0) 



2ai{n • e) + ^ [1 - (n • ef 



(25) 



In Eq. (25) and below we ignore terms of order O (fo/w^) . 
The first term in Eq. ( p5| ) is purely stochastic, because 
(n • e) = 0. Since ((n • e)^) — \/d, we may write this 
stochastic term in the form ri{n) /Vd- The second term 
in Eq. (^5|) has both deterministic and stochastic compo- 
nents, with the latter being negligible in the long time 
limit. The magnitude of the deterministic component is 
{al/vY{l -1/d) ^ {d- l)v^/Av^d. Thus we obtain the 
Langevin equation. 



dvl 
dn 



(26) 



In one dimension, the deterministic term disappears and 
Eq. ( p6| ) coincides with Eq. (p^. 

Following the same steps as those given after Eq. (|l^), 
we eliminate n in favor of the time to transform the above 
equation to 



d\v 



3/2 



■i{d~l)v^ 1 



dt 



lUd 



13/2 



A^/td 



ri{t). 



(27) 



In this equation, the order of magnitudes of the system- 
atic and stochastic terms on the right-hand side are iden- 
tical. Thus evolves by a biased random-walk pro- 
cess, but one in which the bias and the dispersion are of 
the same scale. This can be seen more clearly by writing 
the underlying Fokker-Planck equation for P{u, i) , where 
u = Following the standard prescription ||l^, this 
Fokker-Planck equation is 



dP 
'dt 



lUd 



di 



3d 



du 



P 



(28) 



' — V At h + 



(22) 



with respect to the center of the transparency sphere. 
Here h and e are the unit vectors in the direction of mo- 
tion after the scattering event and the electric field, re- 
spectively, i. e., Vn — vh and E — Ee. The next collision 
event takes place on the surface of the sphere = i'^. 
Consequently, the time increment r between collisions is 
implicitly given by 



(wr)^ -f- avT^{n ■ e) -|- 



aV4 



(23) 



We apply scaling to solve this equation. Let us make the 
scaling ansatz 



P{u,t) 



1 



-f{z) with z = u/{u). 



(29) 



Substituting Eq. ( p^ ) into the Fokker-Planck equation 
(p8|), writing the time and velocity derivatives in terms 
of the scaling variable, the partial differential equation 
can be separated into two ordinary differential equations. 
For the time dependence of (u), we obtain 



The velocity change between collisions is found from 
energy conservation 



dt 



9vl 

sed' 



(30) 
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This then gives a characteristic speed which is propor- 
tional to (uQt/£)^/^ or {a'^ity/^. For the dependence on 
the scahng variable, the scaling function obeys the ordi- 
nary differential equation 



/(z)-z/'(z) = /"(z) + 



d-1 



f{z) f'{z) 



, (31) 



where the prime denotes differentiation with respect to 
z. One integration gives 



3z 



-z]f{z) + A, 



(32) 



where A is a constant. Since f{z) and its first derivative 
vanish faster than any power of z for z oo, A = 0. 
The solution to the resulting equation is 



2(4-d)/6 

r((rf + 2)/6) 



^{d~l)/3 ^-z'/2^ 



(33) 



The numerical coeflScient is determined by the normaliza- 
tion condition f{z)dz = 1, r(y) denotes the gamma 
function |2l|. Notice that the existence of the scaling 
form for the SDF hinged on the magnitude of the bias 
vanishing as u~^. 



D. Distributed Mean-free Paths 

In both the random walk approach for d = I and the ef- 
fective medium theory for d > 1, a mean- free path which 
has the fixed value £ for each scattering event was an 
inherent feature. However, in the Boltzmann equation 
approach of Piasecki and Wajnryb a Poisson distri- 
bution of mean-free paths is implicitly assumed. More- 
over, there will be a distribution of mean-free paths in 
any real scattering medium. We therefore consider the 
physical effects that such a distribution has on trans- 
port properties. We consider a power law distribution of 
mean-free paths, 



p{£) - A'^-i/r, 



(34) 



since this form, for /j, = 3, corresponds [p2|-p4[ to the 
Lorentz gas in a scattering medium with an "infinite" 
horizon (e. g., a square lattice of small scatterers). Prob- 
ability theory ||25| , ^6| suggests that if the distribution is 
relatively sharp, the previous random walk arguments 
apply, while for a broad distribution, different transport 
behavior arises. 

Let us first consider the effect of distributed mean-free 
paths in one dimension. That is, a new mean-free path is 
independently chosen from the above distribution after 
each scattering event. We allow /i to be arbitrary, since 
this general situation is tractable. If the second moment 
of p(£) is finite, i. e., /i > 3, then the distribution of a 
sum of a large number of independent random variables, 



each distributed according to p{£), approaches a Gaus- 
sian and the random walk argument of Sec. II applies. 
In contrast, for ii < 3, a Levy distribution, whose index 
depends on p ||2^,|2^, emerges. Making use of well-known 
results ^5 26 1 for Levy distributions, we determine the n- 
dependence of v^^s to be (the analog of Eq. (||) , 



y/ n Inn; /i = 3, 



(35) 



with Vq = Xa. Repeating the calculational steps in Sec. 
II, we find, for the time dependence of Uims(0 



s(i) 



vo X 



^ Vpt ^ 2,1-3 . 



2 < p<3. 



(36) 



The average displacement in the field direction is thus 
given by {x{t)) = Wi-,„s(t)^/2a, while the drift velocity 
is Vdiih{t) = {x{t))/t. For p, < 2, the first moment of 
p{£) diverges, so that the typical mean-free path is infi- 
nite. Consequently, collisions become irrelevant asymp- 
totically, so that the typical velocity should grow linearly 
in time and an asymmetric velocity distribution should 



IV. NUMERICAL SIMULATIONS 

To test our theoretical predictions, we perform Monte- 
Carlo simulations of particle motion in a two dimensional 
effective medium. An important element in this simula- 
tion is determining where an arbitrary parabolic trajec- 
tory, which starts at the origin, intersects the circumfer- 
ence of a concentric circle. This involves the unwieldy 
solution of a quartic equation. However, since the typi- 
cal speed grows with time, individual trajectory segments 
should deviate only slightly from straight lines, especially 
in the long time limit. Thus we compute the trajectory 
and the time between collisions in a perturbation series 
appropriate for the large speed limit. If the speed hap- 
pens to fall below a preset threshold such that a strongly 
curved trajectory segment should arise, we impose the 
constraint that for this segment the particle is deflected 
exactly parallel to the field. Since the typical speed in- 
creases with time, this "reflecting" boundary condition in 
velocity space is anticipated to have a negligible influence 
on the long-time motion of the test particle. 

Our simulation algorithm therefore consists of the fol- 
lowing basic steps to compute the velocity and time in- 
crements between collisions. These steps are repeated to 
generate a single particle trajectory: 

• If the speed is above a predetermined threshold, 
vth, then: 

1 . Choose a random scattering angle in the range 
0<9 <TT (see Fig. 1). 
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2. Determine the angular position a of the parti- 
cle when it hits the surface of the circle. From 
elementary mechanics, the angle a is pertur- 
batively given in the large velocity limit by 

a = e-esme + e^sm2e + ..., (37) 

with e = {vo/2v)^. 

3. From the angle a, determine the change in the 
longitudinal position of the particle, Ax, and 
thereby determine the change in the speed of 
the particle by Vj = vf + v'^Ax/i. Here Vi is 
the speed of the particle as it begins from the 
center of the transparency circle, and Vf is the 
particle speed just before the collision at the 
circumference of the circle. 

4. Determine the time increment, t, associated 
with this trajectory. In the large velocity 
limit, T is perturbatively given by 

T=-< 1 -ecos6i+ — (cos^6l- 1) + ... L 

(38) 

Here r and v have been expressed in units of 
i/vo and vq respectively. 

• If the speed is less than Vth, then, the scattering an- 
gle is taken to be ^ = 0. Consequently, v'j — vf+vl 
and T = 2. 

Clearly, the sharp difference in the update of the par- 
ticle motion for initial speed smaller or larger than the 
threshold is a crude approximation. One can straightfor- 
wardly construct more accurate, but more cumbersome 
rules to integrate over low-speed trajectory segments. 
However, since these segments are relatively unlikely, this 
refinement was not pursued. Because of the arbitrari- 
ness in the integration of the low-speed segments, the ac- 
tual value of ?;th is also somewhat arbitrary and we chose 
''^th — vq. This appears to provide a relatively good com- 
promise between accuracy and limiting the range over 
which the arbitrary reflecting boundary condition is im- 
posed. To appreciate the numerical implications, con- 
sider, for example, v = Vth- For this case, the maximum 
possible deviation between the initial and final angular 
position of the trajectory arises when 9 « 111.5°, with 
a — 9 ^ —15.8°. For v = 2vth, the maximum deviation 
point occurs when 9 « 97°, with a — 9 ^ —3.7°. Thus 
the effect of the curvature in the individual trajectory 
segments should typically be small. 

We also performed a more faithful simulation for two 
dimensions in which the correct hard-circle scattering is 
implemented. In place of step 1 given above, we as- 
sume that just before the n**^ collision, with incidence 
angle a„_i (see Fig. 1), the test particle uniformly illu- 
minates the cross-section of the scatterer which is taken 



to be a circle of radius r. After specular reflection by 
the scatterer, the difference between the incident and fi- 
nal angles is dip — tt ~ 2sin^^(6/r), where the impact 
parameter b is uniformly distributed between ±r. This 
angular deflection is used to compute the outgoing an- 
gle 9n = (Xn-i + dip, and the corresponding incoming 
angle a„. Our simulation results for this more faith- 
ful implementation of hard-circle scattering are virtually 
identical to those from isotropic scattering. Because of 
this agreement, and also for simplicity, our simulations 
concentrated on the case of isotropic scattering. 
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FIG. 2. Monte Carlo simulation results for 2000 walks of 
1.5^® steps in a two-dimensional effective medium. Shown are 
Wrms(t) (A), and the mean longitudinal position {x{t)) (x). 
The straight lines represents the best fits to the data in the 
range 1.5" < t < 1.5^^ 



Typical results from this Monte Carlo simulation with 
isotropic scattering are presented in Fig. 2. Shown are 

Vn-nsit) and {x(t)) on a double logarithmic scale based on 
2000 trajectories of 1.5^^ « 127,834 steps for the case 
where the transparency circle has a fixed radius. Af- 
ter some transient behavior, the data for t > 500 ap- 
pear to be linear and a linear least-squares fits yields the 
respective slopes of 0.329 and 0.665, in excellent agree- 
ment with the respective theoretical predictions of 1/3 
and 2/3. For true hard-circle scattering in two dimen- 
sions, our simulations gave the corresponding exponent 
estimates of 0.330 and 0.662. Thus the effect of the 
antipersistence in the particle motion truly appears to 
be irrelevant. In Fig. 3, we present corresponding re- 
sults for the distribution of u = at t = 1.5^° (o), 
t = 1.523 (o), t = 1.5^6 (V), and t = 1.529 {+). The raw 
data has been scaled so that the abscissa is z = u/{u), 
while the ordinate is f{z) = {u)P{u,t). This scaled 
data at later times has then been smoothed by aver- 
aging over a small neighborhood to reduce fiuctuations. 
These data compare well with the theoretical prediction 
f{z) = (2V3/r(2/3)) X ^i/3e-^V2 (pig. 3). 
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FIG. 3. Scaled distribution {u)P{u,t) versus the variable 

z = u/{u), where u = |«|^^'^. Representative data shown in- 
clude t = 1.5^" {o),t = 1.5^^ {<>), smoothed over a 3-site neigh- 
borhood, t — 1.5^^' (V), smoothed over a 5-site neighborhood, 
and t = 1.5^® smoothed over a 7-site neighborhood. The 
curve is the theoretical prediction 0.930 . 



V2 



As discussed previously, a lattice array of scatterers 
leads to a power-lavi^ distribution of mean-free paths. 
We therefore also performed simulations of the effective 
medium where the radius of the next transparency cir- 
cle is chosen from the distribution p{£) ^ A^^^/£^, with 

= 3. In this case, we found that the time dependence 
of Wrms(i) and {x{t)) is quite close to that obtained for 
the case of a fixed-radius transparency circle. The data 
for t > 1000 appear to be linear on a double logarithmic 
scale, and a linear least-squares fits in this range yields 
the respective slopes of 0.340 and 0.671 (Fig. 4). For 
this case, however, the data for Wrms(i) and {x{t)) exhibit 
a slight downward trend, a feature which could be at- 
tributed to a logarithmic correction. However, our data 
are insufficient to test for such a correction quantitatively. 
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FIG. 4. Monte Carlo simulation results for 2000 walks of 
1.5^^ steps in a two-dimensional effective medium in which 

the radius £ of the transparency sphere is drawn from the dis- 
tribution p{£) oc £^'^. Shown are Wrms(i) (A), and the mean 
longitudinal position {x{t)} (x). The straight lines represents 
the best fits to the data in the range 1.5^'' <t< 1.5^^. 

The distribution of speeds also exhibits relatively good 
data collapse, but there are quantitative discrepancies 
between the shape of the scaling function and the pre- 
diction /(z) « 0.930 ... X z^/-"^ e""^'/2 ^^la^ ^^^^ ja^^g^ ^j. 
the case of a fixed-radius transparency sphere (Fig. 5) 




FIG. 5. Scaled distribution {u)P{u,t) versus the variable 
z = u/{u), where u = The radius £ of the transparency 

sphere is drawn from the distribution p{£) oc £~^ . Representa- 
tive data shown include t = 1.5^° (o), t = 1.5^^ (o), smoothed 
over a 3-sito neighborhood, t = 1.5^^ (V), smoothed over a 
5-site neighborhood, and t = 1.5^^ (-I-), smoothed over a 7-site 
neighborhood. The curve is 0.930 . . . x z^/^e"^^/^. 



V. DISCUSSION AND SUMMARY 

We have investigated the transport of a charged parti- 
cle which is driven by a constant field in a dissipationless 
elastic and isotropic scattering medium - the field-driven 
Lorentz gas. A fundamental and intriguing feature of this 
system is that the transport is non-stationary. Namely, 
as a function of time, the typical velocity grows as t^/^. 
Although this growth is unbounded, it is significantly 
slower than a linear time dependence that would occur 
in the absence of scattering. In one dimension, we have 
developed a random walk approach in which the parti- 
cle hops to the right or left with equal probability at 
each scattering event. The time increment associated 
with each hop between neighboring scatterers is position 
and direction dependent, a feature which underlies the 
anomalous time dependence of the typical velocity and 
mean displacement. 

Based on this random walk picture, we also obtained 
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the velocity distribution by first writing the Langevin 
equation for the evolution of the typical velocity and the 
underlying Fokker-Planck for the velocity distribution. 
We also constructed a Lifshitz tail argument which re- 
produced the correct behavior for the velocity distribu- 
tion. The solution to the Fokker-Planck equation yields 
the Gaussian distribution in the variable u — \v\^/^, 

P(u,i)(x -i=e-"'/*, (39) 

which, when written in terms of \v\, becomes 

P(i;,i)oc ^He-I^l'/*. (40) 

Interestingly, this is similar, but not coincident with the 
asymptotic velocity distribution function 

P(t;,i)(x^e-HV*. (41) 

obtained by Piasecki and Wajnryb iQ from the Boltz- 
mann equation approach. However, this approach im- 
plicitly assumes an "annealed" medium in which there 
is a Poisson distribution of distances between collision 
events. Thus, while our random walk and the Boltz- 
mann approach are expected to give the same scaling of 
the typical speed with time, the form of the velocity dis- 
tribution from the two approach should not be expected 
to coincide. 

Our random walk argument can also be applied to the 
interesting case of an alternating electric field E{t) = 
Eq sin^ut) to give the counterintuitive result that the 
combination of an AC field and isotropic scattering leads 
to unbounded growth in the speed. This growth arises 
precisely because of the isotropy in the scattering events. 
When the time between collisions becomes less than the 
time for the field to reverse, then the direction of the 
field becomes irrelevant. Consequently, our random walk 
argument for a DC field directly applies and v^^s grows 
without bound in time. The validity of this statement 
depends only on the existence of a well-defined typical 
magnitude for the field. Thus for an AC external elec- 
tric field, the scatterers assist in the absorption of field 
energy by the test particle. In contrast, in an AC field 
with no scattering, a test particle merely follows the field 
and the typical speed is bounded. 

In higher dimensions, we introduced an effective 
medium approximation which provides a simple and 
physically appealing description for the motion of a 
charged test particle. This approximation posits that 
the test particle moves a fixed radial distance along a 
parabolic field-biased trajectory within a "transparency" 
sphere and that an isotropic collision event occurs when 
the particle reaches the surface of this sphere. The as- 
sumption of scattering when a particle moves a fixed 
radial distance implies an annealed medium. Thus one 
might expect a relatively close connection between the ef- 
fective medium and the Boltzmann equation approaches. 



However, because of the bias in the free-particle trajec- 
tory segments, an initially isotropic distribution of out- 
going particle directions immediately after one scattering 
event becomes anisotropic at the next scattering. Within 
an equivalent one-dimensional random walk description 
of the test particle motion, this anisotropy can be de- 
scribed in terms of an effective bias which is proportional 
to The logical consequences of this feature again 

leads to a typical speed which again grows as t^^"^, just 
as in one dimension. 

The effect of the field-induced bias is more pronounced 
in the behavior of the speed distribution, however. Fol- 
lowing a similar approach as that given for one dimen- 
sion, the solution to the Fokker-Planck equation for 
u — \v\^/^ is 

, / \ (d-l)/3 

which, when written in terms of \v\ gives 

^(-'*)- (43) 

while the corresponding result of Piasecki and Wajnryb 
is 

P(.,i)cx^e-H^/*. (44) 

While these two forms agree for d = 2, the coincidence is 
unexpected. The Boltzmann approach explicitly builds 
in isotropy in the collision events and in the interven- 
ing particle motion, while the effective medium explic- 
itly accounts for the field-induced bias between scattering 
events. 

An attractive aspect of the effective medium approach 
is that it can be easily generalized to a distribution of 
mean-free paths, a feature which arises in a lattice re- 
alization of the Lorentz gas. Such a distribution may 
be accounted for by a power-law distribution of sphere 
radii p{£) oc with ^ = 3. This represents a marginal 
case between the regime where distributed radii appear 
to have no effect, for /i > 3, to the case where the scaling 
of the mean speed with time is affected, for 2 < /i < 3. 
Our numerical simulations indicate that the case of /i = 3 
leads to behavior similar to that of no dispersion in the 
sphere radii. However, the applicability of either the 
Boltzmann equation approach or our effective medium 
description to a lattice realization of the Lorentz gas has 
yet to be tested. 

We thank R. S. Chivukula for a helpful discussion 
and J. Machta for particularly useful advice about hard- 
sphere scattering and related suggestions. This research 
was supported in part by the NSF (grants DMR-9219845 
& DMR-9632059), and by the ARO (grant DAAH04-93- 
G-0021). 
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